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SYSTEM FOR COMPUTING 
PROBABILITY DISTRIBUTION OF LOAN LOSSES 




BACKGROUND OF THE INVENTION 
5 Field of the Invention 

This invention relates to a computer technology for 
computing the probability distribution of loan losses which may 
emanate from the entire loans in a financing organization having 
a plurality of loan debtors on condition that the loaned amount 
10 to each loan customer and the bankruptcy probability of each 
customer are known. 

Description of the Related Background Art 

Each of financing organizations of larger scales , such as 
city banks, has as many as tens of thousand loan customers it 

15 finances. Focusing at one loan customer, there are expected two 
cases, namely, it will go into bankruptcy and it will not go into 
bankruptcy. That is, two possibilities of bankruptcy and 
non- bankruptcy can be expected for each loan customer . Therefore , 
when the probability distribution of loan losses emanating from 

20 the entire loaned sum is computed, the entire number of patterns 
or aspects of bankruptcy and non-bankruptcy will amount to a number 
as enormous as self multiplication of 2 over the times corresponding 
to the niomber of customers , that is , the tens of thousandth power 
of 2. Therefore, even when bankruptcy probabilities indicating 

2 5 probabilities of bankruptcy of individual loan customers and loan 

amounts to individual loan customers are already known, 
calculation must be done for all of such enormous patterns which 
are the tens of thousandth power of 2 (approximately one thousandth 
power of one trillion), and computation of probability 
30 distribution is actually impossible even when using a computer. 
That is, even by using computers, their computation load is too 
heavy to actually execute computation. Currently, therefore, 
there is no technology for precisely computing probability 
distributions of loan losses in this field. 

3 5 Since no technology for accurately computing probability 

distributions of loan losses existed heretofore as explained above , 
a simulation method was employed to obtain probability 
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distributions of loan losses. This simulation method was the 
technique of forecasting the probability distribution of the 
entire loan losses by extracting a part of the whole bankruptcy 
and non-banJcruptcy, and calculating probability distributions 
5 of loan losses of these samples. That is, this was the method 
for selecting about ten thousands at random from patterns amounting 
the tens of thousandth power of 2, calculating an aspect of the 
probability distribution of loan losses of these selected patterns , 
and presuming the whole probability distribution of loan losses 

10 involving all of the patterns amounting to the tens of thousandth 
power of 2. A graph of a whole probability distribution of the 
loan losses obtained by this simulation method is shown in Fig. 
12 . As apparent from the graph of Fig . 12 , this simulation method 
was not suitable for practical use because of large errors . That 

15 is, there was a problem that a probability distribution of one 
loan loss could not be obtained because of the low accuracy of 
the graph. In other words, there was a problem that it was not 
suitable for practical use because estimation errors were too 
large. 

20 

SUMMARY OF THE INVENTION 
Taking those problems into consideration, it is an object 
of the invention to provide a computer technique for precisely 
computing probability distribution of loan losses which may occur 

2 5 from the entire loans when loan amounts to individual loan customers 

and bankruptcy probabilities of these loan customers eire known. 
That is, it is the object of the invention to enable precisely 
obtaining probability distribution of loan losses of all loan 
customers or a group of loan customers of a f inancing-^organization 
30 and to enable precisely obtaining a probability density to which 
loan losses are likely to occur. 

According to an aspect of the invention , there is provided 
a system for computing probability distribution of loan losses 
in a financing organization having a plurality of loan customers , 

3 5 comprising: 

an input means for inputting loan amounts loaned to each 
loan customer and bankruptcy probabilities of each loan customer; 
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a characteristic function calculating means for 
calculating a characteristic function on the basis of each loan 
amount and each bankruptcy probability inputted through the input 
means ; 

5 a probability distribution calculating means for 

calculating a probability distribution by Fourier transform 
inversion of the characteristic function calculated by the 
characteristic function calculating means; and 

a probability distribution output means for outputting the 
10 probability distribution calculated by the probability 
distribution calculating means. 

According to another aspect of the invention, there is 
provided a system for computing probability distribution of loan 
losses in a financing organization having a plurality of loan 
1 5 customers , comprising : 

an acquiring means for acquiring rating fluctuation 
probabilities which are the probabilities of fluctuation in 
ratings of each loan customer, and credit value changing amounts 
which are changes in value of the credit to each loan customer 
20 caused by the fluctuation in the ratings thereof; 

a characteristic function calculating means for 
calculating a characteristic function on the basis of each rating 
fluctuation probability and each credit value changing amount; 

a probability distribution calculating means for 
25 calculating a probability distribution by Fourier transform 
inversion of the characteristic function calculated by the 
characteristic function calculating means; and 

a probability distribution output means for outputting the 
probability distribution calculated by the probability 
30 distribution calculating means. 

According to a further aspect of the invention, there is 
provided a system for computing probability distribution of loan 
losses in a financing organization having a plurality of loan 
customers , comprising : 
35 an input means for inputting loan amounts loaned to each 

loan customer and bankruptcy probabilities of each loan customer; 

a actual loss calculating means for calculating actual 
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losses which are amounts of the losses the financing organization 
may actually suffer when each loan customer goes into bankxuptcy; 

a characteristic function calculating means for 
calculating a characteristic function on the basis of the actual 
5 losses and the bankruptcy probabilities; 

a probability distribution calculating means for 
calculating a probability distribution by Fourier transform 
inversion of the characteristic function calculated by the 
characteristic function calculating means; and 
10 a probability distribution output means for output ting the 

probability distribution calculated by the probability 
distribution calculating means. 

According to a still further aspect of the invention, there 
is provided a system for computing probability distribution of 
15 loan losses in a financing organization having a plurality of 
loan customers, comprising: 

a scenario acquiring means for acquiring a plurality of 
actual losses which are amounts of losses the financing 
organization may suffer when each loan customer goes into 
20 bankruptcy by predicting future fluctuations, and a plxirality 
of bankruptcy probabilities of each loan customer by predicting 
future fluctuations, and using these values as a plurality of 
scenarios ; 

a characteristic function calculating means for 
2 5 calculating characteristic function for each scenario on the basis 
of the actual losses and the bankruptcy probabilities acquired 
by the scenario acquiring means; 

a probability distribution calculating means for 
calculating probability distributions for each scenario by Fourier 
30 transform inversion of the characteristic function calculated 
by the characteristic function calculating means; 

an average probability distribution calculating means for 
calculating an average probability distribution which is the 
average of the probability distributions for each scenario; and 
35 a probability distribution output means for outputting the 

average probability distribution calculated by the average 
probability distribution calculating means . 
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According to a yet further aspect of the invention , there 
is provided a system for computing probability distribution of 
loan losses in a financing organization having a plurality of 
loan customers, comprising: 
5 a scenario acquiring means for acquiring a plurality of 

actual losses which are amounts of losses the financing 
organization may suffer when each loan customer goes into 
bankruptcy by predicting future fluctuations, and a plurality 
of banlcruptcy probabilities of each loan customer by predicting 
10 future fluctuations, and using these values as a plurality of 
scenarios ; 

a characteristic function calculating means for 
calculating characteristic functions for each scenario on the 
basis of the plurality of actual losses and bankruptcy 
15 probabilities acquired by the scenario acquiring means; 

a probability distribution calculating means for 
calculating probability distributions for each scenario by Fourier 
transform inversion of the characteristic function calculated 
by the characteristic function calculating means; and 
20 a probability distribution output means for outputting the 

probability distributions calculated by the probability 
distribution calculating means for each the scenario. 

According to a further aspect of the invention, there is 
provided a system for computing probability distribution of loan 

2 5 losses in a financing organization having a plurality of loan 

customers , comprising : 

a scenario acquiring means for acquiring loan amounts of 
each loan customer and acquiring a plxjrality of bankruptcy 
probabilities of each loan customer by predicting future 
30 fluctuations, and then using these values as a plurality of 
scenarios ; 

a characteristic function calculating means for 
calculating characteristic functions for each scenario on the 
basis of the loan amounts and the bankruptcy probabilities acquired 

3 5 by the scenario acquiring means; 

a probability distribution calculating means for 
calculating probability distributions for each scenario by Fourier 
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transform inversion of the characteristic function calculated 
by the characteristic function calculating means; 

an average probability distribution calculating means for 
calculating an average probability distribution which is the 
5 average of the probability distributions for each scenario; and 

a probability distribution output means for outputting the 
average probability distribution calculated by the average 
probability distribution calculating means. 

According to a still further aspect of the invention, there 
10 is provided a system for computing probability distribution of 
loan losses in a financing organization having a plurality of 
loan customers, comprising: 

a scenario acquiring means for acquiring loan amounts of 
each loan customer and acquiring a plxirality of banlcruptcy 
15 probabilities of each loan customer by predicting future 
fluctuations, and then using these values as a plurality of 
scenarios ; 

a characteristic function calculating means for 
calculating characteristic functions for each scenario on the 
2 0 basis of the loan amounts and the bankruptcy probabilities acquired 
by the scenario acquiring means; 

a probability distribution calculating means for 
calculating probability distributions for each scenario by Fourier 
transform inversion of the characteristic function calculated 
2 5 by the characteristic function calculating means; and 

a probability distribution output means for outputting the 
probability distributions calculated by the probability 
distribution calculating means for each the scenario. 

According to a yet further aspect of the invention, there 
30 is provided a system for computing probability distribution of 
loan losses in a financing organization having N loan customers 
k=l.,.N comprising: 

an input means for inputting loan amounts Mj^ to individual 
N loan customers k=l . . .N and bankruptcy probabilities Pj^ thereof ; 
35 a loan customer calculating means for calculating the number 

N of the loan customers on the basis of the loan amounts 1%. and/or 
the bankruptcy probabilities Pj^ inputted by the input means; 



a characteristic fianction calculating means for 
calculating a characteristic function 

N 



0(O = JJ{t + i^.(exp(rWj-l)} 



5 at each t of t=2 7rm/(2^") (m=0, 1, 2 2^"-l) for the number 

of points n of Fourier transform; 

a probability distribution calculating means for 
calculating a probability distribution by Fourier transform 
inversion of the characteristic function calculated by the 
10 characteristic function calculating means , by using a fast Fourier 
transform technique; and 

a probability distribution output means for outputting the 
probability distribution calculated by the probability 
distribution calculating means. 
15 According to a yet further aspect of the invention, there 

is provided a method for computing probability distribution of 
loan losses in a financing organization having a plurality of 
loan customers, comprising: 

an input process for inputting loan amounts loaned to each 
20 loan customer and bankruptcy probabilities of each loan customer; 

a characteristic function calculating process for 
calculating a characteristic function on the basis of the loan 
losses and the bankruptcy probabilities; 

a probability distribution calculating process for 
25 calculating a probability distribution by Fourier transform 
inversion of the characteristic function calculated by the 
characteristic function calculating process; and 

a probability distribution output process for outputting 
the probability distribution calculated by the probability 
30 distribution calculating process. 

According to a yet further aspect of the invention, there 
is provided a recording medium storing a program for calculating 
a probability distribution of loan losses in a financing 
organization having a plurality of loan customers, in which the 
35 program comprising: 

an input process for inputting loan amounts loaned to each 
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loan customer and bankruptcy probabilities of each loan customer; 

a characteristic function calculating process for 
calculating a characteristic function on the basis of the loan 
losses and the bankruptcy probabilities; 
5 a probability distribution calculating process for 

calculating a probability distribution by Fourier transform 
inversion of the characteristic function calculated by the 
characteristic function calculating process; and 

a probability distribution output process for outputting 
10 the probability distribution calculated by the probability 
distribution calculating process. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a diagram showing the relation between Fourier 
15 transform and Fourier transform inversion of a probability 
distribution and a characteristic function; 

Fig. 2 is a flow chart generally showing a process for 
obtaining a probability distribution of loan losses according 
to the first embodiment of the invention; 
2 0 Fig. 3 is a flow chart showing the process for obtaining 

the probability distribution of loan losses according to the first 
embodiment in greater detail ( 1 ) ; 

Fig. 4 is a flow chart showing the process for obtaining 
the probability distribution of loan losses according to the first 

2 5 embodiment in greater detail (2); 

Fig. 5 is a diagram showing a graph of a probability 
distribution obtained by the invention; 

Fig. 6 is a flow chart generally showing a process for 
obtaining a probability distribution of loan losses according 
30 to the second embodiment of the invention; 

Fig. 7 is a diagram showing a graph of a probability 
distribution obtained by the invention; 

Fig. 8 is a diagram for explaining patterns of possible 
changes in rating of a certain customer from the current time 

3 5 up to a certain future time; 

Fig, 9 is a flow chart generally showing a process for 
obtaining a probability distribution of loan losses according 
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to the third embodiment of the invention; 

Fig. 10 is a flow chart generally showing a process for 
obtaining a probability distribution of loan losses according 
to the foirrth embodiment of the invention; 
5 Fig . 11 is a system for computing a probability distribution 

of loan losses, which realizes the invention by hardware; and 

Fig. 12 is a diagram showing a graph of a probability 
distribution obtained a conventional simulation method. 

10 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

(First Embodiment) 

The first embodiment of the invention is configured to first 
obtain a characteristic function of loan losses of the whole loan 
customers in a financing organization, relying on the feature 

15 that a probability distribution and a characteristic function 
are related by one-to-one in Fourier transform and Fourier 
transform inversion, and thereafter obtain the probability 
distribution of the whole loan losses by Fourier transform 
inversion of the cheiracteristic function . This is explained below 

20 in greater detail. 

First referring to Fig. 1, explanation is made about the 
relation between Fourier transform and Fourier transform inversion 
related to the probability distribution and the characteristic 
function. Fig. 1 is a diagram showing a relation between Fourier 

25 transform and Fourier transform inversion of a probability 
distribution and a characteristic function. 

As shown in Fig, 1, by Fourier transform of the probability 
distribution, the characteristic function is determined. In 
contrast, by Fourier transform inversion of the characteristic 

30 function, the probability distribution is determined. The 
embodiment shown here uses this relationship to first obtain a 
characteristic function of the entire loan losses and then obtain 
the probability distribution of the entire loan losses by Fourier 
transform inversion of the characteristic function obtained. 

35 Next explained is a process for obtaining the probability 

distribution of the entire loan losses according to this 
embodiment . 
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Let the number of loan customers be N. That is, the number 
of loan customers for obtaining the probability distribution of 
loan losses is N. The number of loan customers N may be either 
the number of all loan customers of a certain financing organization 
5 or the number of all loan customers within a group of collected 
plural loan customers. Xj, is the random variable of bankruptcy 
(Xj,=l) or non- bankruptcy {X^=0) of a loan customer k. Here, k=l, 
2, . . . N. The bankruptcy probability which is the probability 
that the loan customer k goes into bankruptcy is p,^. Further, 

10 let the random variable be independent for each loan customer. 
That is, it is assiamed that the probability that a certain loan 
customer goes into bankxuptcy and the probability that another 
loan customer goes into bankruptcy are independent from each other . 
When the loan amount to the loan customer k is Mj, and the loan 

15 loss in the entire loan amount is L, distribution of the loan 
loss L in the entire loan amount can be expressed by Equation 
(1). 



N 



(1) 



20 



25 



Consider the characteristic function 0 (t) of the loan loss 
L in the entire loan. Definition of the characteristic function 

( t ) is characteristic function (/) ( t ) =E Cexp ( iLt ) ] where i= V-T , 
E is the expected value, and exp() is the exponential function 
having the base of a natural logarithm as its base. The 
characteristic function (f>(t) can be developed as follows. 



<p(t) = E[txp(iLt)] 



exp 



30 



Because of the independence. 
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N 



E[txp(itMkXk)] 



N 



{(1 - pk) + pk cxp(itMk)} 




(2) 



From Equation (2), the characteristic function can be 
formulated. 

Next consider a Fourier transform. Here is used fast 
Fourier transform ( FFT ) using a computer . By this FFT , when complex 
nimbers are given for the values of the characteristic function 
<f)(t) in 



then, Fourier transform inversion takes place, and the probability 
distribution f ( L ) of the loan loss L in the entire loan is computed . 
Here, m=0, 1, 2, . . . 2"-l, andL=0, 1,2, ... 2"-l. n is the parameter 
for determining the number of points Q for Fourier transform. 
That is, the number of points Q = 2^". Therefore, by sequentially 
substituting Equation ( 3 ) into Equation ( 2 ) , the probability 
distribution f (L) of the loan loss L can be computed. 

Next referring to Figs . 2 through 4 , a process for obtaining 
the probability distribution f (L) of the loan loss L mentioned 
above is explained. Fig. 2 is a diagram generally showing the 
process flow for obtaining the probability distribution f (L) of 
the loan loss L. Figs. 3 and 4 are diagrams showing the process 
flow in greater detail. 

First explained with reference to Fig, 2 is an outline of 
the process executing the technique for computing the probability 
distribution of a loan loss according to the embodiment through 
a computer . As shown in Fig . 2 , first inputted are the loan amounts 
to each loan customer and their probabilities Pj, of bankoruptcy 
(step SI). Subsequently, after computing the number of loan 
customers N by using at least one data of these loan amounts Mj, 
and bankruptcy probabilities Pj,, the characteristic function <f) 
(t) is computed (step 82) . Then, by conducting Fourier transform 
inversion of this characteristic function (f>(t) , the probability 




(3) 
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distribution f (L) of the loan loss L in the entire loan amount 
is obtained ( step S3 ) , And, the probability distribution is output 
through a printer, for example (step S4). This is the outline 
of the process by this embodiment. 

Next made with reference to Fig . 3 is a detailed explanation 
about the process executing the technique for computing the 
probability distribution of the loan loss according to the 
embodiment through a computer • As shown in Fig . 3 , the loan amounts 
to each loan customer and their probabilities p,, of bankruptcy 
are inputted as data ( step Sll ) . k is the variable for identifying 
the loan customer. That is, M(k) is the loan amount to the loan 
customer k, and P{k) is the bankruptcy probability of the loan 
customer k. If these data are downloaded from a large-scale host 
computer of a bank, for example, they can be inputted more easily. 
Next, as apparent from Equation (3) , while fixing t=0 (step S12) , 
1 is selected for the variable k for identifying the loan customer, 
and 1 is selected for the variable A for obtaining the characteristic 
function (t> {t) corresponding to a certain value of t (step S13) . 

Next conducted an arithmetic operation based on Equation 
(2) shown above. That is, A=AX {i+p(k) X (exp(t XM(k) ) -1} is 
repeated sequentially from k=l for the number of loan customers 
N (step S14 through step S16). Subsequently, the computed A is 

substituted for the characteristic function 0(t) . Then, 2 7r/(2^) 
is added to t (step S18). That is, the next t is computed. If 
this t does not reach 271 (step S19) , the process is repeated from 
step S13. Vftien t reaches or goes beyond 27C (step S19) , it means 
that the characteristic function 0 (t) has been obtained, and the 
process of Fourier transform inversion shown in Fig. 4 is executed. 

That is, as shown in Fig. 4, a Fourier coefficient is 
extracted from the characteristic function <!> (t) (step S20), 
Subsequently, by dividing 2^" as the number of points Q by 4, the 
number of points Q is quartered (step S21). Then, Fourier 
coefficients of 4Xq points are regarded as Q sets of serial 
four-point Fourier coefficients (step S22). That is, regarding 
a series of coefficients C(w) , if the remainder of division w 
by Q is w', C(w' ) , C(w'+Q), C(w'+2Q) and C(w'+3Q) are regarded 
as a series of four-point Fourier coefficients. After that. 
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four-point Fourier transform inversion is repeated Q times (step 
S23 ) . That is , four-point Fourier transform inversion is executed 
for Q sets . These series of coefficients obtained are expressed 
as 

5 D(S, w') S=0, 1, 2 and 3 

Thereafter, exp( ( -2 7riw'S/(4Q) ) which is a rotation factor 
is multiplied on Q sets of four-point Fourier coefficients series 
{stepS24) . Then, Q sets of four-point Foiirier coefficients series 
calculated in this manner are regarded as four sets of Q-point 
10 Fourier coefficients series (step S25) , and it is judged whether 
Q=l is satisfied or not (step S26). If it is judged in step S26 
that Q=l is not satisfied, Q is divided by 4 (step S27) , and the 
process is repeated from step S21. If it is judged in step S26 
that Q=l is satisfied, it means that the Fourier transform inversion 

15 of the characteristic function (t> it) has been completed and the 
probability distribution has been obtained. Therefore, putting 
the probability density in the probability distribution on the 
ordinate and putting the loan loss on the abscissa, the probability 
distribution is outputted (step S28) . An example of this output 

20 of the probability distribution is shown in Fig. 5. The graph 
shown in Fig. 5 is made by drawing loan losses on the abscissa 
in bar graphs in the unit of one yen. 

As understood from Fig. 5, this embodJLment calculates the 
probability distribution only up to loan losses beyond which the 

25 probability density of loan losses can be regarded substantially 
zero in the calculation process . That is , the probability density 
of loan losses above a predetermined amount may be regarded zero 
either in the calculation process or in the practical use. 
Therefore, by omitting calculation of the probability distribution 

30 of loan losses above a predetermined amount , the embodiment intends 
to increase the speed of computation. For example, in Fig. 5, 
the probability density of loan losses above four hundred billion 
yen can be regarded substantially zero, this embodiment obtains 
the probability distribution only up to loan losses of four hundred 

35 billion yen. In this case, n=20. That is, since a scale for four 
hundred billion yen in the unit of one yen is required as the 
abscissa of Fig. 5, n may be selected so that 2^''-l exceeds at 
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least four hxindred billion. That is, since 2^^^°-l is about one 
trillion and nine hundred ninety five billion, n=20 may be set. 

It is also possible to calculate probability distribution 
about all loan losses without omitting those above a predetermined 
5 amount unlike this embodiment. That is, calculation process may 
be done so as to output graphs up to the total amount of loans 
to all loan customers. 

According to the technique of calculating the probability 
distribution of loan losses according to the embodiment, the 

10 probability distribution of loan losses can be computed accurately 
by using a computer when there are a number of loan customers. 
Therefore, financing organizations like banks can previously 
estimate occurrence of loan losses as a probability . For example , 
taking the example of Fig . 5 , it can be expected as a probability 

15 that the probability density of the loan loss of 70 billion yen 
is 0.0010. Also, it can be known that the probability density 
of loan losses above three hundred billion yen is 99%. By 
previously knowing occurrence of loan losses as a probability, 
banks and other financing organization can do risk management 

20 on credits by preparing reserve for loan loss or completing their 
own capitals. 

Additionally, since probability distribution of loan 
losses of and above which the probability density can be regarded 
substantially zero is omitted from calculation, the load to a 
2 5 computer can reduced , and probability distribution can be computed 
in a short time. 



(Second Embodiment) 

The second embodiment of the invention determines a 
30 predetermined amount like one million yen or one hundred million 
yen as a unit, then rates respective loan amounts of individual 
loan customers into groups based on the predetermined unit, and 
calculating probability distributions for respective groups to 
alleviate the computation load to the computer. 
3 5 Fig. 6 is a diagram showing a process of computing 

probability distribution of loan losses in a financing 
organization according to the second embodiment in form of a flow. 
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As shown in Fig. 6, the process flow of the second embodiment 
is different from the process flow of the first embodiment shown 
in Fig. 2 in that step S30 is interposed between step SI and step 
S2 . That is , after inputting the loan amounts and the bankruptcy 
5 probabilities p,^ (step SI) concerning each loan customer, these 
loan amounts are rounded to the place of the predetermined unit , 
that is, rounded into integer times of the predetermined unit 
{ step S30 ) . In this rounding process , loan amoxints J\ are adjusted 
by rounding up, rounding down or rounding off to the place of 

10 the predetermined unit, i.e. , to numbers which are integer times 
of the predetermined unit. In this embodiment, loan amounts 
are rounded up to numbers which are integer times of one hundred 
million yen , After that , the value of n necessary for computation 
is determined. That is, an amount of loan loss above which the 

15 probability of occurrence of loan losses can be regarded 
substantially zero is divided by the predetermined unit, one 
hundred million yen , to obtain its quotient . Then , n is determined 
so that the maximum value of m in Equation (3) , namely, 2^-1 be 
equal to or larger than that quotient. Similarly to the first 

20 embodiment, assuming that the amount above which the probability 
of loan losses can be regarded substantially zero is four hundred 
billion yen, for example, if the predetermined minimum unit is 
ten million yen, its quotient is 4000. Therefore, the number of 
points Q for Fourier transform becomes 4000, and since 2^^-1=4095, 

25 n is determined as 6. 

It is also possible to calculate probability distribution 
about all loan losses instead of limiting it to a certain amount 
of loan loss, unlike this embodiment. In this case, by dividing 
the total amount of loans by a predetermined minimum unit, and 

30 n may be determined so that 2^"-l be equal to or larger than the 
quotient . 

The process after step S30 is the same as that of the first 
embodiment. In this embodiment, the graph output of the 
probability distribution of loan losses is bar graphs of every 
35 levels corresponding to integer times the predetermined minimum 
value. An example of the graph output is shown in Fig. 7. As 
shown in Fig. 7, the graph output shows bar graphs on the abscissa 
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graduated in one hundredyen increments of the loan losses . However , 
by using a smaller unit amount, an output with substantially no 
error can be obtained. 

According to the technique of calculating the probability 
5 distribution of loan losses according to this embodiment, the 
time required for calculation can be shortened by alleviating 
the load to the computer. That is, although the first embodiment 
can calculate probability distribution over the range of loans 
from 0 to 2^"-!, it has the problem that it uses one yen as the 

10 minimum unit and n becomes unacceptably large in a financing 
organization having the total amount of loans which reaches or 
exceeds one trillion. The second embodiment, however, can 
decrease the value of n by rounding loans of loan customers to 
niambers which are integer times of a predetermined unit. As a 

15 result, this embodiment can alleviate the calculation load to 
the computer and shorten the calculation time. Therefore, this 
embodiment can be realized by using an inexpensive personal 
computer, for example, and the cost for developing a system can 
be reduced - 

20 Additionally, in this embodiment, when loans to loan 

customers are rounded up to numbers which are integer times of 
a predetermined unit, they are rounded up. Therefore, a graph 
of probability distribution taking the worst estimation into 
account can be obtained. That is, a graph of probability 

2 5 distribution with a fail-safe concept can be obtained. 

(Third Embodiment) 

The third embodiment of the invention is directed to 
calculation of probability distribution of broadly defined loan 

30 losses, considering changes in value of credits caused by changes 
in the "rating", which shows credible levels of loan customers 
(debtors) , as a kind of loan losses in a broader sense. That is, 
although the first and second embodiments take only two states, 
i.e. whether loan customers go or do not go into bankruptcy, to 

35 consider that loan losses occur only upon bankruptcy, the 
embodiment shown here considers that loan losses will occur when 
the rating of loan customers changes even when they do not go 



17 



into bankruptcy. 

With reference to Fig, 8, explanation is made about the 
"rating". In the example of Fig. 8, there are seven grades of 
rating of loan customers , "AAA", "AA" , "A", "BBB", "BE", "B" and 
"D" . "AAA" is the highest grade , and indicates that a loan customer 
graded so is a blue-chip enterprise. The rating sec[uentially 
degrades to "AA" , "A", "BBB" , "BB" , and "B", and "D" indicates 
that a loan customer so graded went into bankruptcy. 

A change in value of credits caused by a change in "rating" 
means that, if all credits are re-evaluated at a reference point 
of time in a future, credits to loan customers changed to higher 
ranks will be evaluated higher and credits to loan customers changed 
to lower ranks will be evaluated lower . This embodiment considers 
changes in value of credits to be involved in loan losses in a 
broad sense, and intends to comprehend in form of probability 
distribution how the value of credits is expected to change at 
a future point of time. 

In the example of Fig. 8, a specific loan customer is 
currently ranked "A" . The probability P,^ that this loan customer 
k is ranked "AAA" at a future point of time is 1%, the probability 
Pj^ of being ranked "AA" is 5%, the probability Pj^ of being ranked 
"A" is 80%, the probability P,^ of being ranked "BBB" is 6%, the 
probability Pj^ of being ranked "BB" is 5%, the probability P^ 
of being ranked "B" is 2%, and the probability P,^ of being ranked 
"D" is 1%, Note here that, also when the loan customer k remains 
ranked "A" at the future point of time, it is expressed in 
probability P^, assuming for convenience that it changes from 
the current grade "A" to the grade "A" . 

Here is expected that, if the loan customer is ranked "AAA" 
at the future point of time, the ratio of the loss of the credit 
value caused by the change of the rating will be -10% of the loan 
amount. That is, it is expected that a profit of 10% of the loan 
amount will be generated. If the loan customer is ranked "AA" 
at the future point of time, the loss ratio of the credit value 
is expected to be -5% of the loan amount. That is, a profit of 
5% of the loan amount is expected. If the loan customer is ranked 
"A" at the future point of time, there is no change in the rating. 
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and the credit value is expected to remain unchanged. 

If the loan customer k is ranked "BBB" at the future point 
of time, the loss ratio of the credit value caused by the change 
in the rating is expected to be 5% of the loan amount. That is, 
a loss of 5% of the loan amount is expected. If the loan customer 
k is ranked "BE" at the future point of time, the loss ratio of 
the credit value caused by the change in the rating is expected 
to be 10% of the loan amount. That is, a loss of 10% of the loan 
amount is expected. If the loan customer is ranked "B" at the 
future point of time, the loss ratio of the credit value caused 
by the change in the rating is expected to be 20% of the loan 
amount. That is, a loss of 20% of the loan amount is expected. 
If the loan customer k is ranked "D" at the future point of time, 
it means bankruptcy of the loan customer, and the loss of all 
loan amount is expected . 

Next, the characteristic function (t> (t) taking such a 
possible change in rating into account is derived. 

Assume that there are N loan customers k=l, 2, . . . N, and 
H grades of rating. Let Pj^ be the probability that the rating 
of the loan customer k changes to an h-order grade. Since the 
loan customer k has one of grades of rating. 



Let Mj^ be the changed amount of the credit value as a result 



of a change of the rating of the loan customer k to the h-order 
grade at a certain point of time in the future. In the example 
of Fig. 8, by multiplying the loss ratio by the change in rating 
by the loan amount, the changed amount of the credit value 
can be obtained. Under these conditions, the characteristic 
function 0 ( t ) of the changed amount of the credit value is 





(5) 




19 

That is , by merely modifying Equation { 2 ) to Equation { 5 ) , the 
characteristic function of the changed amount of the credit value, 
i.e. the characteristic fxinction (f>(t) of the loan losses in a 
broad sense can be expressed. 
5 By processing the characteristic function ( t ) by Fourier 

transform inversion, probability distribution of loan losses in 
a broader sense can be obtained. 

Next referring to Fig. 9, explanation is made about a 
technique for calculating probability distribution of the 

10 broad- sense loan losses according to this embodiment through a 
computer. Fig. 9 is a flow chart for explaining the process 
therefor. As shown in Fig. 9, loan amounts M^, probabilities p,^ 
of a change in rating and loss ratios caused by the change in 
rating are inputted first (step S40) . The probabilities p,^ of a 

15 chcinge in rating and the loss ratios caused by the change in rating 
may be inputted for every loan customer or they may collectively 
determined for all possible changes in rating independently from 
specific loan customers . 

Next, after computing the number of loan customers N by 

20 using the data of loans M^, the changed amount in credit value 
for each loan customer is calculated from the loan of each loan 
customer and the loss ratio caused by a change in rating (step 
S41). 

After that , using the probability P,^ of the change in rating 
2 5 of each loan customer and the changed amount of the credit 

value of each loan customer, the characteristic function 0(t) 
shown in Equation (5) is computed (step S42). After that, by 

Fourier transform inversion of the characteristic function <f) 
(t), the probability distribution f(L) of the broad- sense loan 

30 losses is obtained ( step S3 ) . Then, this probability distribution 
is output through a printer, for example (step S4). 

As explained above, the computing technique concerning 
probability distribution of broad- sense loan losses according 
to the invention enables risk management of credits taking account 

35 of the risk that each credit value may decrease due to a change 
in credit the loan customer toward the worse even when it does 
not go into bankruptcy. That is, as shown in Fig. 8, even if a 
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loan customer does not fall into bankruptcy (grade "D" ) at a future 
point of time, if its rating falls to the grade "BBB", "BE" or 
"B" , the value of the credit thereto decreases. This decrease 
in value of the credit can be regarded to be a loan loss in a 
broad sense. The third embodiment can obtain probability 
distribution of broad- sense loan losses, taking such decrease 
of the credit value into consideration. 

(Fourth Embodiment) 

The fourth embodiment of the invention is directed to 
precisely obtaining probability distribution of loan losses by 
preparing a plurality of scenarios about actual losses and 
probability of loan losses and averaging these scenarios. 

The actual loss herein means the amount obtained by 
subtracting a secured amount , etc, from a loan to a loan customer. 
That is, it can be expressed by actual loss = (loan to the loan 
customer) - (secirred amount) . This is because, even when a loan 
customer goes into bankruptcy, if a guarantee or collateral is 
held for the loan, a part or all of the loan as much as the value 
of the collateral can be recovered, and an amount corresponding 
to the difference is the actual loss. 

First explained is a general formula. Assume that there 
are N loan customers k=l, 2, . . . N, their bankruptcy probability 
is Pj,, and probability that the loan amount is ^ is G (p^, . . . , 
Pn, . . . , . And assume that probability distribution of loan 
losses obtained by Fourier transform inversion is f (L, p^, . . . , 
Pj^, . . . , M^) in the case where the bankruptcy probability of 

each loan customer is p^, . . . , p^,, and the loan amount is M^, . . . , 
. In such cases where the bankruptcy probability and loan amounts 
vary with a certain probability, probability distribution of loan 
losses becomes 

^ f(L, pi, pN, Ml, Mn)G(pi, pN, Ml, Mn) ( 6 ) 





35 



(where siams are calculated on all combinations of p^, p^, M^, 
'•\ M^) 
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Alternatively, in the case where probability distribution 
of bankruptcy probability and changes in loan amounts is continuous , 
it is assumed that there are N loan customers k=l, 2, . . . N, their 
5 bankruptcy probability is p,^, and the probability density is 
g(Pj^, . . . , Pj^, . . . , ^^^) . In this case, probability distribution 
becomes 



J*f(I^ pi, pN, Ml, MN)g(pi, pN, Ml,..., MN)dpi...dpNdMi...dMN ( 7 ) 

10 

That is, by numerical integration, probability distribution of 
loan losses is obtained. 

Next referring to Fig. 10, explanation is made about a 
process executing the technique for calculating probability 
15 distribution of loan losses according to this embodiment by means 
of a computer. As shown in Fig. 10, loan amounts ^ to individual 
loan customers and their bankruptcy probabilities p^ (step SI) 
are inputted (step S50) . Next , after computing the number of loein 
customers N by using at least one of the data of loans sn scenarios 
20 are generated, and the actual loss of each loan customer and its 
bankruptcy probability are generated for each scenario . That is , 
about each of the first to sn-th scenarios, the actual loss of 
each loan customer and its bankruptcy probability are generated. 

For example, if the value of the collateral changes in the 

2 5 future, it results in a change of the actual loss. The bankruptcy 

probability may also change with possible economical fluctuation 
in the future. As to actual losses and bankruptcy probabilities 
for these sn scenarios , operator ' s expected values maybe inputted, 
or they may be calculated by computer simulation or by using a 
30 function. 

After that, probability distribution of loan losses is 
calculated for each scenario. That is, the first actual losses 
and the first bankruptcy probabilities of individual loan 
customers according to the first scenario are taken in ( step S52 ( 1 ) ) , 

3 5 and the characteristic function is computed based on these first 

actual losses and first bankruptcy probabilities (step S53(l) ) . 
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Subsequently , the characteristic function are processed by Fourier 
transform inversion ( step S54 { 1 ) ) to obtain the first probability 
distribution (step S55{1)). 

Similarly, the second probability distribution according 
to the second scenario is acquired (stepS52(2) through stepS55( 2) ) , 
the third probability distribution according to the third scenario 
is acquired (step S52{3) through step S55(3) ) , et. seq. , and the 
sn-th probability distribution according to the sn-th scenario 
is acquired (step S52(sn) through step S55(sn)). 

Through these steps, sn probability distributions are 
obtained. In this embodiment, based on Equation (6) shown above, 
a weighted mean of these sn probability distributions is obtained 
(step S56) , and the probability distribution as the weighted mean 
is output through a printer, for example (step S47). Although 
the example of Fig. 10 obtains the probability distribution of 
loan losses averaging all scenarios, it is also possible to 
individually output sn probability distributions of loan losses 
through a printer, for example, so that an operator can analyze 
the outputted results . Additionally, if Equation (7) shown above 
is used, average of the scenarios is obtained by numerical 
integration in step 856. 

Next explained is a technique for generating a plurality 
of scenarios based on a function in step S51. In this embodiment, 
a plurality of bankruptcy probabilities of respective loan 
customers are generated on the basis of a function. Therefore, 
regarding the loan amounts, they may be fixed for respective loan 
customers , or data predicted by an operator for respective 
scenarios may be input. 

(Multi-factor Model) 

First explained is a multi-factor model. Let Nonn( ) be 
a cumulative distribution function of standard normal 
distribution . Assume that there are R random variables u^ , U2 , ... 
Ur , they are in accordance with an R- dimensional normal distribution , 
and its probability density function is fR( ) . Also assume that 
the number of loan customers is N, that there are random variables 
y,^ which indicate conditions of respective loan customers k=l , . . . , 
N, and that the random variables can be expressed by 
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(8) 



where a,^ are constants, are random variables independent from 
the random variables Ui, Uj, . . . Up in accordance with the standard 
normal distribution, and respective 6 are independent from each 
other. Further assume that there are constants Y,,, and when yk<Y„, 
a loan customer goes into bankxuptcy. 

In this case, the possibility that the loan customer k goes 
into bankruptcy, that is, 

R 

are equal in meaning, and bankruptcy and non- bankruptcy of the 
loan customer k is determined only by the random variable £ under 
the condition where the random variables u^, U2, . . . Ur are fixed. 

Since random variables s are independent from each other, 
bankruptcy and non -bankruptcy of the loan customer k are 
independent from each other under the condition where the random 
vciriables u^, Uj, . . . u^ are fixed. In this case, the bankruptcy 
probability of the loan customer k is specified as 



Norm 



(9) 



That is, from Equation (9), probability distribution of 
R- dimensional bankruptcy probabilities can be obtained . From the 
R-dimensional function expressed by Equation (9), a necessary 
number of bankruptcy probabilities are sampled. That is, sn 
scenarios are acquired from the R-dimensional function expressed 
by Equation ( 9 ) , and bankruptcy probabilities of individual loan 
customers are acquired for each of the scenarios . More 
specifically, dmn values are sampled about each factor u^. of 



r 

y » 
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Equation (9), and by repeating it for R dimensions, {dinn)^=sn 
scenarios are acquired. Then, for each scenario, bankruptcy 
probabilities of individual loan customers are acquired from 
Equation ( 9 ) . 

5 In this example of this embodiment, sn is 25 ♦ The random 

variable u^ is a factor indicating the general economical condition , 
and the smaller the value, the worse the economical condition. 
R pieces are prepared as the factor indicating the economical 
condition to predict various economical conditions, a,^ is a 
10 coefficient indicating in the r dimension how the business of 
the loan customer k is liable to suffer influences from economical 
fluctuations, and the larger the value, the more liable to suffer 
influences from the general economical condition the loan customer 
k. 

15 If the probability distribution of loan losses is F (L, 

u^, U2, ... Ur) , theaveragedprobability distribution of loan losses 
is 



jF(L,u^,'",Uj^)f^{u^,'",Uj^)du^-"duj^ (10) 

20 

By numerically integrating it, the probability distribution of 
loan losses can be obtained. 

(One -factor Model) 
25 As a special case of the multi-factor model, there is a 

one-factor model . Next made below is explanation on the one-factor 
model . 

LetNorm( ) be a cumulative distribution function of standard 
normal distributions. Assume that a random variable u exists, 
30 and it is accordance with a standard normal distribution. Also 
assume that the number of loan customers is N, there are random 
variables y,, which indicate conditions of respective loan customers 
k=l, N, and the random variables can be expressed by 
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y„=a,,u+ £ J, 



(11) 
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where a^^ is a constant, random variable independent from 

the random variable u in accordance with the standard normal 

distribution, and respective £^ are independent from each other. 
Further assume that there is a constant Yj,., and when yk,<Yi^# a loan 
5 customer goes into bankruptcy. 

In this case, the possibility that the loan customer k goes 
into bankruptcy, that is, 

yK<YK 

10 £K<Yk-a,,u 

are equal in meaning, and bankruptcy and non- bankruptcy of the 
loan customer k is determined only by the random variable s under 
the condition where the random variable u is fixed - 

15 Since random variables £ are independent from each other, 

bankruptcy and non-bankruptcy of the loan customer k are 
independent from each other vinder the condition where the random 
variable u is fixed. In this case, the bankruptcy probability 
of the loan customer k is specified as 



20 



Nortn{y,,-a„u) (12) 



That is, from the function expressed by Equation (12), sn 
probabilities of bankruptcy are sampled. That is, sn scenarios 

25 are acquired from the function expressed by Equation (12), and 
bankruptcy probabilities of individual loan customers are acquired 
for each of the scenarios. 

In this embodiment, sn is 25. The random variable u is 
a factor indicating the general economical condition, and the 

30 smaller the value, the worse the economical condition, a,^ is a 
coefficient indicating how the business of the loan customer k 
is liable to suffer influences from economical fluctuations, and 
the larger the value, the more liable to suffer influences from 
the general economical condition the loan customer k. 

35 If the random distribution of loan losses is F (L, u) , the 

averaged probability distribution of loan losses is 
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Sf(L, u)Nonn'(u)du 



(13) 



By numerically integrating it, the averaged probability 
distribution of loan losses can be obtained. 

As explained above , the technique for computing probability 
distribution of loan losses according to the embodiment can obtain 
more precise probability distribution of loan losses by preparing 
a plurality of scenarios different in actual loss and bankruptcy 
probability, then calculating a plurality of probability 
distributions based on these scenarios and obtaining their average . 
That is, it is possible to obtain probability distribution of 
loan losses taking account of possible future general economical 
aspects and fluctuations in value of the collateral. 



probabilities and fluctuations in loan amounts, it is useful for 
risk management of credits in a financing organization. 

The invention is not limited to the above -explained 
embodiments, but can be modified and changed in various modes. 
For example, probability density of loan losses can be output 
in form of a numerical list instead of a bar graph. Further, 
although the foregoing embodiments have been explained as using 
fast Fourier transform for executing Fourier transform inversion , 
other techniques, such as the use of typical Fourier transform 
inversion formulas, are also acceptable. That is, by selecting 
Q as the number of points, Fourier transform inversion can be 
executed by using the Fourier inverse transform formula shown 
as Eqpiation (14). 



value of a collateral or guarantee held for a loan as explained 
in the fourth embodiment , probability distribution of loan losses 
may be computed by using the actual loss which a financing 
organization will actually lose when the loan customer goes into 



Since this technique can calculate probability 
distribution of loan losses upon fluctuations in bankruptcy 




(14) 



In the first to third embodiments , taking account of the 
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bankruptcy . 

Furthermore, it is also possible to record the processes 
according to the above -explained embodiments on a recording medium 
such as floppy disk or GD-ROM. In this case, by retrieving the 
5 data from the storage of the recording medium onto a general -purpose 
computer , the technique according to the invention can be realized . 

It is also possible to realize those embodiments in a 
hardware mode . Configuration of an exemplary system for computing 
probability distribution of loan losses realizing the first 

10 embodiment, for example, in a hardware mode is shown in Fig. 11. 
As shown in Fig. 11, the system for computing probability 
distribution of loan losses includes a loan amount and bankruptcy 
probability input device 10 through which loan amounts and 
bankiruptcy probabilities of individual loan customers are inputted, 

15 a characteristic function computing device 12 for computing a 
characteristic function on the basis of these loan amounts and 
bankruptcy probabilities, a probability distribution computing 
device 14 for computing probability distribution by Fourier 
transform inversion of the characteristic function, and a 

20 probability distribution output device 16 for output ting the 
calculated probability distribution as a graph to a printer. 

According to the present invention, a financing 
organization having a plurality of loan customers can precisely 
compute probability distribution of loan losses through a computer 

25 and stochastically estimate occurrence of loan losses in advance- 



